Solid-solid collapse transition in a two dimensional model molecular system 
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Solid-solid collapse transition in open framework structures is ubiquitous in nature. The real dif- 
ficulty in understanding detailed microscopic aspects of such transitions in molecular systems arises 
from the interplay between different energy and length scales involved in molecular systems, often 
mediated through a solvent. In this work we employ Monte Carlo (MC) simulations to study the 
collapse transition in a model molecular system interacting via both isotropic as well as anisotropic 
interactions having different length and energy scales. The model we use is known as Mercedes-Benz 
(MB) which for a specific set of parameters sustains three solid phases: honeycomb, oblique and tri- 
angular. In order to study the temperature induced collapse transition, we start with a metastable 
honeycomb solid and induce transition by heating. High density oblique solid so formed has two 
characteristic length scales corresponding to isotropic and anisotropic parts of interaction poten- 
tial. Contrary to the common believe and classical nucleation theory, interestingly, we find linear 
strip-like nucleating clusters having significantly different order and average coordination number 
than the bulk stable phase. In the early stage of growth, the cluster grows as linear strip followed 
by branched and ring-like strips. The geometry of growing cluster is a consequence of the delicate 
balance between two types of interactions which enables the dominance of stabilizing energy over 
the destabilizing surface energy. The nuclei of stable oblique phase are wetted by intermediate order 
particles which minimizes the surface free energy. We observe different pathways for pressure and 
temperature induced transitions. 



I. INTRODUCTION 

Several studies have shown that many solids exhibit 
temperature and pressure induced solid-solid phase tran- 
sitions resulting in a discontinuous change in material 
properties [il-Hl- In the case of shock wave induced solid- 
solid transitions both temperature and pressure of the 
system increases [IjQ- The well known examples of solid- 
solid transition are Fe and its alloys (Fe/Ni alloys trans- 
form on cooling from face-centered cubic (fee) phase to 
a low-temperature body-centered cubic (bcc) phase), al- 
loys based on CuAl, NiTi, NiAl, and oxide ceramics (like 
Zr02), name only a few. In steel, the structural tran- 
sitions are used to enhance its strength. Carbon crys- 
tallizes in the hexagonal graphite structure at ambient 
pressure and undergoes a pressure-induced transition to 
the tetrahedral diamond structure. Structural transfor- 
mations are also observed in biological systems. Some 
virus species use the pressure-induced martensitic trans- 
formation to infect bacteria cells [1|. Another applica- 
tion of solid-solid transition is the shape-memory effect 
in many alloys and widely used in medicinal science 3]. 
Solid-solid transitions can be driven either by homoge- 
nous or heterogeneous (in the case of presence of preex- 
isting defects which reduces the free energy barrier for 
transition) nucleation. Despite the great importance of 
solid-solid transitions in material, medicinal and biologi- 
cal sciences, little is known about the microscopic aspects 
of solid-solid transitions. In molecular systems, the real 
difficulty in understanding detailed microscopic aspects 
of such transitions arises from the interplay between dif- 
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ferent energy and length scales involved and often medi- 
ated through a solvent. 

Computer simulations are the most appropriate tech- 
nique to study the microscopic aspects of solid-solid tran- 
sitions. However, the difficulty in computer simulations 
arises in construction of potential which can give rise to 
multiple solid phases. Now it is well accepted that sys- 
tems interacting with core soft potentials show rich phase 
diagram consisting multiple solid phases |8l-[l3|. Core- 
softened potential was first used by Hemmer and Stell^ 
in a lattice gas system to discuss the isostructural solid- 
solid phase transition. Later, Young and Alder 9;] showed 
that a two dimensional system interacting via repulsive 
step potential (hard-core plus square-shoulder repulsion) 
can rise to phase diagram and solid-solid phase transi- 
tions qualitatively similar to those observed in cesium 
and cerium. Distinct solid phases in systems interacting 
with core soft isotropic potential arise due to the pres- 
ence of distinct length scales in the potential. Recently, 
Sengupta and coworker |13J introduced a potential where 
particles interact via modified Lennard- Jones potential 
and each particle has two internal states. They also ob- 
served a rich phase diagram showing water-like thermo- 
dynamic anomalies and polymorphism. Although there 
are many computer simulation studies that discuss the 
rich phase diagram consisting many solid phases |8l-[l4| 
only few studies are devoted for understanding of nature 
and pathway of solid-solid transitions. Recently, Rao and 
Sengupta[15|-[l7| studied the nucleation of a triangular 
phase from bulk metastable square crystal in a model two 
dimensional system interacting via both two and three 
body potentials. These studies suggest that at high tem- 
perature the critical nucleus of the triangular phase is 
isotropic and at low temperature the critical nucleus is 
anisotropic twinned bicrystal. They also observed that 



the final microstructures depend on quench rate. For sys- 
tems like iron the dependence of microstructures on rate 
of cooling (quench depth) is well known. Slow cooling 
rate results in an equilibrium ferrite phase while rapid 
cooling results the trapping of system in a metastable 
twinned martensite. 

In this work we have studied the microscopic pathway 
for the solid-solid collapse transition in a model two di- 
mensional (2D) molecular system interacting via both 
isotropic as well as anisotropic interactions having differ- 
ent length and energy scales. The model we use is known 
as Mercedes-Benz (MB) potential |18l - l20{ which on spe- 
cific parameterization for relative strength of isotropic in- 
teraction sustains three solid phases: honeycomb, oblique 
and triangular [2l!|. The complexity of MB model and 
simplicity of 2D allows us to explore several aspects of 
solid-solid transitions not observed in atomic systems. In 
the subsequent sections we shall discuss the honeycomb 
to oblique solid collapse transition. 

The organization of the rest of the paper is as follows. 
In the next section (Section II), we present the model 
system and the simulation details used in this work. In 
Sections III, first, we discuss the free energy surface and 
then the microscopic pathway for temperature and pres- 
sure induced honeycomb to oblique solid transition. In 
Section IV we present concluding remarks. 



which is a measure of relative strength of the isotropic 
to anisotropic interaction is defined as A,; = £lj/£hb- 
In the present study we have used \i — 0.4. Note that 
we have two length as well as energy scales correspond- 
ing to the isotropic LJ and anisotropic HB interactions. 
Monte Carlo simulations are performed in constant pres- 
sure and temperature (NPT) ensemble having number 
of particles N = 576, 836 and 1008. All quantities are 
represented in reduced units scaled with HB parameters 
such as temperature as ksT /\£hb\ and distance is scaled 
with hydrogen-bond length, uhb ■ 

In the case of temperature induced transition (except Fig. 
IC) we started with honeycomb structure at T = 0.05 
and then increased the temperature in steps of 0.01 up 
to T = 0.10 and then increased in steps of 0.0005 at P 
= 0.10. We first equilibrated the system for 5 x 10^ MC 
steps and then collected the trajectories for next 5 x 10^ 
MC steps. For pressure induced transition we started 
with honeycomb lattice at pressure P = 0.10 and then 
increased the pressure in steps of 0.05 up to P = 0.4 and 
then increased in steps of 0.01 at T = 0.05. In this case 
the equilibration step was 1.5 x 10^ MC steps and then 
collected the trajectories for next 1 x 10^ MC steps. 



III. RESULTS AND DISCUSSION 



II. MODEL AND SIMULATION DETAILS 

Mercedes-Benz potential has two terms - isotropic 
Lennard- Jones (LJ) and an anisotropic hydrogen bond 
(HB) term written in a completely separable form as [isl - 



U{X,,X.j) = ULj{nj) + Uhb{X.,Xj), 



(1) 



where Xi denotes vectors representing both the coordi- 
nates and orientation of the i*^ particle and rij the dis- 
tance between the centers of two molecules. Ulj is de- 
fined as 

UlMj) = ^eLj[{<JL.j/njf^ - {cTL.j/n.j)% (2) 

where Elj and ulj are the wel-depth and diameter 
for the isotropic LJ interaction. The anisotropic HB 
part,t//fB, is defined as 

UHB{Xi,Xj) = eHBG{rij — ghb) 

■^ (3) 

X ^ G{ik.u,, - l)G{ji.u^j + 1), 

k,l=l 

where G{x) = exp {-x'^ /2(t'^) with a = 0.085. The unit 
vector ik represents the fc*'' arm of the i*'' particle (fc — 
1, 2, 3) and Uij is the unit vector joining the center of 
particle i to the center of particle j. The HB parame- 
ters are Ehb = — 1 and bond length aHB = 1. LJ con- 
tact distance {(Tlj) is O-lanB- Anisotropic parameter 



A. Free energy surface 

As the form of the Mercedes-Benz potential is writ- 
ten as a completely separable isotropic and anisotropic 
parts having distinct ground states (triangular and hon- 
eycomb), one would expect that the interplay between 
the relative strength of the isotropic to anisotropic inter- 
action can give rise distinct solid phases - either stable 
or metastable. Recently, we find that the interplay be- 
tween the strength of isotropic and anisotropic interac- 
tions gives rise to rich phase diagram[21]. The computed 
phase diagram consists of isotropic liquid and three crys- 
talline phases - honeycomb, oblique and triangular. In- 
terestingly, we note that the oblique phase is neither the 
ground state of isotropic nor the anisotropic part of the 
potential. It forms due to delicate balance between the 
two types of interactions. In order to study the solid- 
solid collapse transition from low density open structured 
honeycomb to high density oblique phase the honeycomb 
solid (ground state of the anisotropic part of the poten- 
tial) should be metastable (as shown in Fig. lA) at 
a given temperature, pressure and relative strength of 
isotropic interaction (Xi). We have constructed such a 
free energy surface by tuning the relative strength of the 
isotropic interaction (Xi). 

The honeycomb lattice is stabilized by 3 hydrogen 
bonds (HBs) per particle and the oblique phase is sta- 
bilized by two HBs per particle along with LJ interac- 
tions (see Fig. IB). Alignment of bonds in Fig. IB indi- 
cates hydrogen bonding. Since both the phases are solid 
and stabilized by anisotropic interactions (either com- 
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FIG. 1. (A) The proposed free energy surface for honeycomb to oblique solid-solid transition at Xi = 0.4 is shown. (B) The 
structures of the honeycomb and oblique solid phases are shown. Alignment of bonds indicates hydrogen bonding. (C) Isotropic 
(ulj) and anisotropic (uhb) interaction energies per particle and the total energy per particle (u) before and after transition 
are shown. Blue line (with triangle up symbol) indicates the dependence of enthalpy change per particle (h) on temperature. 
Note the very weak temperature dependence of energy/enthalpy and switch in the contributions of isotropic and anisotropic 
energies to the total energy on transition. 



pletely or partially), one would expect that entropy of 
two phases will not differ significantly. Only vibrational 
entropy (and some contribution from mobility of defects) 
is present and also not expected to differ significantly as 
both the solid phases are stabilized by anisotropic inter- 
actions, though, the fraction of directional bonding de- 
creases on going from honeycomb to oblique phase. This 
indicates that entropy has very little role (compared to 
gas- liquid and liquid-solid transitions) in the construction 
of free energy surface for solid-solid transitions. Thus the 
temperature change should not alter significantly the free 
energy surface; however, it is important to overcome the 
required free energy barrier for the transition. 
To support the proposed free energy surface (Fig. lA) in 
Fig. IC, we have shown the temperature dependence of 
enthalpy, energy and relative contribution of the isotropic 
(LJ) and anisotropic (HB) interactions to the total en- 
ergy per particle of the system for isotropic to anisotropic 
strength ratio, Xi = 0.4. The starting configuration is 
honeycomb lattice at pressure P = 0.10 and T = 0.05, 



which on heating collapses to the high density oblique 
phases. Although, in Fig. IC we have shown the de- 
pendence of energy and enthalpy on temperature near 
the transition temperature only, we find that very weak 
dependence persists up to much lower temperature T = 
0.05. Note that, contrary to the ordinary solid-liquid 
transition (where total energy of the system increases), 
in this case the total energy of the system decreases. A 
very weak dependence of enthalpy on temperature and 
decrease in the enthalpy of the system on transition sup- 
ports the free energy surface proposed in Fig. lA. Hon- 
eycomb solid is metastable at given initial pressure and 
temperature (T = 0.05 and P = 0.10). We also note that 
the relative contributions of LJ and HB parts of interac- 
tion potential to the total energy per particle is switching. 
This indicates the increased importance of isotropic in- 
teraction in the oblique solid phase. 

In order to gain quantitative insight into the struc- 
tures of honeycomb and oblique phases, in Fig. 2 we 
have shown the computed radial distribution functions 
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FIG. 2. (A) Radial distribution functions for honeycomb and 
oblique solid phases are shown. (B) Bond length distribu- 
tions of nearest neighbors for honeycomb and oblique solid 
phases are shown. Note the splitting of first peak of radial 
distribution function for the oblique and also the bimodal 
bond length distribution of nearest neighbors. The splitting 
of the 1"* peak of radial distribution function and bimodal 
distribution of nearest neighbor bond length distribution sug- 
gest that the oblique solid has two characteristic length scales 
corresponding to hydrogen-bond and vander Waals (LJ) in- 
teractions. 



and bond length distributions of nearest neighbors for 
both honeycomb and obhque solid phases near transition 
temperature. We have used the Voronoi construction in 
order to identify the nearest neighbors (NN) of each par- 
ticle. Note the splitting of first peak of radial distribu- 
tion function and also the bimodal bond length distri- 
bution of nearest neighbors for the oblique solid phase. 
This signifies that in the oblique solid all nearest neigh- 
bor particles are not equidistant. There are two length 
scales present in the system, which is a consequence of 
the two length scale potential. Particles having lower 
bond length are stabilized by isotropic LJ interaction and 
larger bond lengths are stabilized by anisotropic HB in- 
teraction. Contrary to the oblique solid, honeycomb solid 
is stabilized only by anisotropic HB interactions and thus 
the nearest neighbor distance distribution has only one 



peak centered at HB length and there is no splitting in 
the first peak of radial distribution function. 
In the subsequent sections, first, we shall discuss the tem- 
perature induced transition followed by a comparative 
study of pressure and temperature induced transitions. 



B. Order parameters and their susceptibilities 

In Fig. 3, we have shown the order parameters change 
and their susceptibilities for temperature induced col- 
lapse transition. On increasing temperature, thermal 
fluctuations destabilize the honeycomb structure and col- 
lapses to a high density oblique structure with large 
change in density (from ^ 0.79 to ^ 1.56) (see Fig. 3A). 
In the collapsed transition, high density change indicates 
that the transition can only be first order. As due to 
increase in order parameter difference the surface energy 
cost for formation of new phase also increases. Since the 
barrier in the global free energy surface as function of 
global order parameter (density and order) is related to 
the surface tension, the probability of long wavelength 
fluctuations of the global order parameters (density fluc- 
tuation is coupled with order fluctuation) in the system is 
less. In Fig. 3B, we have reported the change of both the 
global six- fold bond orientational order (tl^eg) and local 
three- fold bond orientational order {1JJ31) with tempera- 
ture, m-fold global bond orientational order {ij^mg) is 
defined as 
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and m-fold local bond orientational order (ipmi) is defined 
as 




NN 

NN ^^ 



(5) 



where m = 6 for six-fold and m — 3 for three-fold sym- 
metry. NN is the number of nearest neighbors, 6ij is the 
angle of the bond that the f^ neighbor makes with the 
tagged J*'' particle. Significantly large values of global or- 
der parameter, "i/^gg, in both phases indicate that there is a 
global translation of either six- fold or three- fold bond ori- 
entational order. As tp6g does not distinguish between the 
six-fold and three-fold bond orientation order, we have 
also plotted the local three-fold bond orientational order 
(note that global three-fold bond orientational order is 
zero for both six- fold and three- fold symmetry). How- 
ever, the later does not guarantee the long range trans- 
lation of local bond orientational order. Sharp transition 
in -03/ and significantly large values of ipeg in both the 
phases indicate a one step transition from three-fold hon- 
eycomb lattice to six-fold (distorted) oblique solid phase. 
In Fig. 3C and 3D, we have plotted the susceptibil- 
ities along density (isothermal compressibility, kt) and 
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FIG. 3. (A) The dependence of the density (p) of the system on temperature is shown. The starting configuration is honeycomb 
solid and after transition it goes to the oblique phase. (B) The dependence of global six- fold (V'eg) and local three fold ("03;) 
bond orientational order parameters on temperature is shown. In (C) and (D), we have shown the dependence susceptibilities 
of global order parameters (density and six- fold bond orientational order parameter): isothermal compressibility (kt) and six- 
fold bond orientational order susceptibility (xe) on temperature. Weak temperature dependence and finite jump in response 
functions suggests strongly 1"* order transition. 



global six- fold bond orientational order parameter (xe)- 
Six-fold bond orientational order susceptibility is defined 
as X6 — {\ip6g\^) — (|'/'6g|)^ and and isothermal com- 
pressibility is defined in terms of volume fluctuations as 
KT = i{V^) - {Vf)/T{V). Weak temperature depen- 
dence and finite jump in both kt and xe suggest the 
strongly 1** order transition. 



C. Microscopic pathway of temperature induced 
collapse transition 

In order to find the pathway of transition first we need 
to distinguish the particles belonging to new phase in 
the sea of metastable bulk honeycomb solid phase. To 
distinguish the nucleating and growing new oblique-like 
particles from the parent metastable honeycomb phase, 
we computed the distribution of local 3-fold bond orien- 
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tational order (ipsi), defined as V'si = J2j= 
for both the solid phases. As shown in Fig. 4, the com- 
plete separation of two distributions enables us to sepa- 
rate the growing new phase from the parent honeycomb 
solid phase. Due to significant overlap in the distribu- 
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FIG. 4. Distributions of 3-fold local bond orientational order 
("i/jsi) for honeycomb and oblique solids are shown. 



tion of local six-fold bond orientational order we have 
not used it as a parameter to distinguish the two phases. 
Note that, though, i/j^i is able to separate the oblique 
solid from the honeycomb solid phase, it does not distin- 
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FIG. 5. Snapshots of the system at (A) early (B) intermediate and (C) late stage of transition are shown. (D) Oblique solid 
phase after transition for A^ = 836 particles are shown. Note the linear strip-like arrangement transforming into branched and 
ring-like strips. Oblique solid has string-like arrangement of particles. We also note the wetting of oblique phase by a layer of 
intermediate solid phase. (E) Snapshot of a growing nucleus is shown. 



guish the oblique solid-like and liquid-like particles (as 
the peak of distribution is centered at very low value of 
ip3i). A two order parameter criteria using the values of 
both V'si and ipei is required to separate the oblique phase 
with liquid-like particles as the oblique solid is character- 
ized with low ip3i and high ipei values (see Fig. 3B). 
In Fig. 5 A - 5D, we have shown the snapshots of the 
system at different stages during the collapse transition. 
Depending on ip^i value we have defined three types of 
particles in the system. The particles having -03^ value 
greater than 0.75 is considered as honeycomb type and 
indicated with black color, red color particles have inter- 
mediate order with ip^i value between 0.25 and 0.75. The 
particles having -031 value less than 0.25 and ipsi value 
greater than 0.2 (lower cut-off value of ^gi below which 
the distribution of "061 (-P(V'6i)) in oblique phase is zero) 
are considered as oblique-like. Note that these criteria 



do not distinguish oblique-like and triangular solid-like 
particles. However, simplicity of 2D systems allows us 
to distinguish these two types of orders by visual inspec- 
tion of snapshots. Recently, Onuki et al. developed an 
elegant approach based on local bond orientational order 
parameter to quantify different phases, stacking faults 
and grain boundaries in both two and three dimensional 
systems [H, HI]. 

Contrary to the common believe and classical nucle- 
ation theories that the nuclei of new phase are circular 
in 2D) and have same order and density as of bulk stable 
phase, interestingly, we find highly anisotropic growth 
of new phase. At the early stage of growth, the cluster 
grows as linear strip (having significantly different order 
than the bulk oblique phase) followed by branched and 
ring-like strip and finally engulfs the whole system. The 
geometry of nucleating and growing cluster (shown in 
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FIG. 6. (A) The dependence of local order parameter, l-i/'3!c (black filled circles) and average coordination number scaled with 
six (the average coordination number of particles in oblique solid phase), Nic/G (red filled diamonds) of the largest cluster on 
the size of the largest cluster is shown. For convenience in comparison we have plotted 1 — ^sic in place of f/'sic- AH quantities 
are averaged over 100 MC steps along the growth trajectory for A'' = 836 particle system. We have also observed similar trend 
for N = 1008. The blue dotted line indicates the critical cluster size up to which there is no appreciable change in either 
average order or the coordination number of the cluster. (B) Pathways of solid-solid transition in two order parameter (density 
and six-fold bond orientational order) space are shown. The black filled circles represent the evolution of global six-fold bond 
orientational (tpag) order with density and red filled diamonds represent the evolution of local six- fold bond orientational order 
(ipei) with density. 



Fig. 5E) is a consequence of delicate balance between 
the two (isotropic and anisotropic) types of interactions 
which enables the dominance of stabilizing energy over 
the destabihzing surface energy. Sharing of the H-bonds 
(see Fig. 5E) between the growing cluster and the bulk 
metastable honeycomb phase minimizes the surface en- 
ergy cost for formation of new phase. We also observe 
that the structure of the growing cluster is different from 
the bulk oblique solid phase. As shown in Fig. 5B and 
5C, the oblique-like phase starts appearing inside the 
high density intermediate order cluster. There is no clear 
dividing surface between the oblique and the honeycomb 
and hence it is difficult to assign unambiguously the iden- 
tity of particles as honeycomb or oblique-like. Wetting 
of the oblique-like phase by intermediate order particles 
reduces the surface free energy cost for the formation of 
stable oblique phase. 

Thus, the solid-solid transition in a molecular system is 
a two step process. In the first step small anisotropic 
cluster having significantly different order nucleates and 
grows. Anisotropic shape enables the cluster to maximize 
the energy gain due to strong specific interactions among 
particles inside the cluster as well as with the bulk parent 
honeycomb phase. In the second step, due to internal 
rearrangement inside the dense (intermediate ordered) 
anisotropic cluster oblique-like phases start appearing. 
Wetting of nuclei by intermediate order phase is also ob- 
served in crystal nucleation from metastable liquid [2J| 
where nucleation preferentially takes place in regions of 
high structural order. Note that clear dividing surface 
with no sharing of particles between the two phases (clus- 
ter of stable phase and bulk metastable phase) creates 



a density/order deficiency and thus increases the sur- 
face energy [23. Wetting is one of the factors responsible 
for the observed discrepancy between the predictions of 
classical nucleation theory and experimental/simulation 
studies. 

In Fig. 6A, we have shown the evolution of the aver- 
age order and average coordination number of a growing 
cluster. We consider a particle not as honeycomb like if 
it has V'3j value less than 0.75. Particles having ^3^ value 
less than 0.75 and are connected by neighborhood form 
a cluster of new phase. The order of the cluster {ipsic) is 
defined as, i/'szc = Y^'i=i V'si/'^ic? where nic is the size of 
the largest cluster and summation is over all the particles 
that belong to the largest cluster. Average coordination 
number is scaled with the average coordination number 
of particles in the oblique solid phase. Since only one 
cluster grows in the system the largest cluster contains 
the maximum fraction of particles belonging to the new 
phase. The cluster of new phase has significantly differ- 
ent order from the bulk stable phase. During the growth, 
both the order and average coordination number of the 
cluster gradually evolve and finally merge to the value 
corresponding the bulk oblique solid phase. As shown in 
the figure, there is a crossover (indicated by blue dotted 
line in Fig. 6A) from weakly dependence to non-linear 
dependence of order as well as average coordination num- 
ber on size of the cluster. At early stage of growth, weak 
dependence of the order and average coordination num- 
ber on size of the cluster indicates the linear strip-like 
growth, where surface to volume ratio does not change 
appreciably. The non-linear dependence of order/average 
coordination number on cluster size is a consequence of 
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two factors: (i) the change of surface to bulk ratio due to 
branching and swelhng of the growing cluster and (ii) the 
rearrangement of the particles inside the cluster towards 
oblique-like structure. 

In Fig. 6B, we have shown the evolution of both global 
six- fold bond orientational order {tpeg) and local six- fold 
bond orientational order (ipei) of system with density 
during collapse transition. We observe an almost mono- 
tonic decrease of ^g; during transition from the value in 
the honeycomb to its value in the oblique solid phase. 
However, ipQg shows strongly non-monotonic evolution 
from its value in the honeycomb to the oblique solid 
phase. This non-monotonic evolution of V'eg and mono- 
tonic (or very weakly non-monotonic) evolution of tJjqi 
arises due to cancellation of components of global bond 
orientational order at intermediate stage of growth and 
suggests the presence of domains (separated by grain 
boundaries) having different symmetry axes. Presence 
of domains having distinct symmetry axes does not alter 
the values of local order parameters. This scenario is also 
evident from the snapshots shown in Fig. 5. In the late 
stage of growth due to relaxation of grain boundaries the 
global order parameter value again increases and finally 
merges to the value corresponding to the bulk oblique 
phase. 



D. Does the pathway of pressure induced transition 
differ from temperature induced transition? 

Temperature plays two roles in the phase transition ~ 
(i) it couples with entropy and perturbs the entire free 
energy landscape and (ii) it provides sufficient thermal 
energy to overcome the free energy barrier for phase tran- 
sition. In the case of condensation and crystallization, 
where the entropy difference between the two phases is 
large, the dependence of free energy barrier on temper- 
ature is stronger than the case of solid-solid phase tran- 
sitions. As discussed earlier, in the case of solid-solid 
transition where specific interactions play important role 
in the stability of both the phases we can assume that 
the major role of temperature is to provide sufficient ther- 
mal energy to overcome the free energy barrier for tran- 
sition. As the honeycomb solid is metastable at a given 
initial pressure and temperature with sufficiently high 
free energy barrier, transition can take place in two ways 
- (i) increase the temperature until the thermal energy 
is sufficient to overcome the free energy barrier (this is 
called temperature induced transition and discussed ear- 
lier) and (ii) increase the pressure until the given ini- 
tial thermal energy is sufficient to overcome the barrier. 
Increase of pressure favors high density phase and thus 
decreases the free energy barrier for transition. In this 
section we shall compare these two scenarios and discuss 
how the pathway of pressure induced transition at low 
temperature differ from temperature induced transition. 
In Fig. 7 we have compared the pathways for pressure 
and temperature induced transitions in two dimensional 



order parameter (density and six-fold bond orientational 
order) space. For temperature induced transition the 
transition temperature is T = 0.105 at P = 0.10. In the 
case of pressure induced transition we have increased the 
pressure of the system (details are mentioned in Section 
II) at constant temperature T = 0.05 and the observed 
transition pressure is P '^ 0.45. On increasing pressure 
transition temperature shifts towards lower temperature. 
In the case of pressure induced transition, we find that 
different trajectories collapse to different metastable dis- 
ordered states (with slightly different density and order) 
and transition pressure also vary in orders of 0.01. As 
shown in the figure, contrary to the temperature induced 
transition where the honeycomb solid collapses to high 
density crystalline oblique solid phase, in the case of 
pressure induced transition the honeycomb solid collapses 
to a long lived disordered metastable state having very 
small global order parameter value (see Fig. 7A). Disor- 
dered metastable phase evolves towards the stable state 
with an extremely slow rate. In spite of the absence of 
the global six-fold bond orientational order, significantly 
large value of the local six-fold bond orientational order 
(see Fig. 7B) suggests the heterogeneous structure of the 
metastable phase where significant local six-fold symme- 
try is present but there is no global translation of the 
local order. 

In order to gain more insight in the microscopic struc- 
ture of the disordered solid phase, in Fig. 8, we have 
shown the distributions of local six-fold bond orienta- 
tional order (V'ei) for both the oblique and the disordered 
solid phase. The emergence of the shoulder in the distri- 
bution corresponding to the disordered solid phase indi- 
cates the presence of heterogeneous local environments. 
To find the structures of different local environments re- 
sponsible for the broadening of ipei distribution towards 
the lower value, in inset we have shown the distributions 
of Tpmi (P('0mi)) fo^ ^hc particlcs having coordination 
number m. That is we have computed the distribution 
of Tpei for only those particles which have coordination 
number six and so on. The distribution is normalized 
with the total number of particles so that area under 
curve of P{ipmi) will give the fraction of particles having 
coordination number m. We note the presence of sig- 
nificant five-cordinated (~ 30%) and seven-coordinated 
(~ 18%) along with six-coordinated (^ 50%) particles 
with respective symmetries. On heating the metastable 
disordered solid phase undergoes a transition to the sta- 
ble oblique phase having unimodal i/'ei distribution. On 
increasing pressure the peak of the distribution of ijjQi in 
the oblique solid phase shifts towards lower ipQi value. 
Recently, Kadau et a/.'5i] used molecular dynamics sim- 
ulations to investigate the shock-induced phase transfor- 
mation of solid iron. Above a critical shock strength the 
nucleation and growth become collective as many small 
grains start growing simultaneously. However, we have 
observed a single cluster growth scenario in both pressure 
and temperature induced transitions. Collective nucle- 
ation and growth scenario indicates the absence of free 




FIG. 7. Pathways for pressure and temperature induced transitions in two dimensional order parameter (density and global 
six-fold bond orientational order) space are shown. (A) density (p) and global six-fold bond orientational order parameter 
{4>6g) space. (B) density (p) and local six- fold bond orientational order parameter {4'6i) space. Black filled circles represent the 
pathway for temperature induced transition and red filled diamonds represent the pathway for pressure induced transition. 




FIG. 8. Distributions of local six-fold bond orientational or- 
der (f/'ei) for the collapsed structure in the case of temper- 
ature induced transition (black line distribution with filled 
circles) and pressure induced transition (red line distribution 
with filled diamonds) are shown. In case of temperature in- 
duced transition the collapsed structure is oblique solid and 
in the case of pressure induced transition the collapsed struc- 
ture is disordered solid. In inset, distributions of m-fold local 
bond orientational order {il)mi) of particles having coordina- 
tion number m for the pressure induced collapsed phase are 
shown. In the figure the distribution with black line indicates 
the %l)^i distribution of particles having coordination number 
five, red line indicates ij^ei distribution of particles having co- 
ordination number six and blue line indicates -i/jyi distribution 
of particles having coordination number seven. Area under 
'ijjmi distribution gives the fraction of particles having coordi- 
nation number m. 



energy barrier for nucleation. Note that shock waves in- 
crease both the pressure as weU as temperature and incre- 
ment in both decreases the free energy barrier for tran- 
sition. Above a critical shock strength the free energy 
barrier might be negligibly small and hence many clus- 
ters start growing simultaneously. This transition closely 



resembles with the gas-liquid transition at large metasta- 
bilty (above kinetic spinodal) where one observes similar 
collective nature of nucleation and growth[26 28]. In the 
case of pressure induced transition (or low temperature 
transition as on increasing pressure transition point shifts 
towards the lower temperature) trapping of the system in 
a long lived local minima indicates that the free energy 
landscape is rugged. Although one might observe the 
rugged energy landscape in atomic systems also, inclu- 
sion of specificity in interaction potential enhances the 
ruggedness of the landscape. The origin of enhanced 
ruggedness lies in the competition between specific and 
non-specific interactions where local regions might be en- 
ergetically stabilized by a particular type of interaction 
but system as a whole might not be in global minimum. 
These local traps in the landscape decreases significantly 
the dynamics of the transition. On increasing tempera- 
ture the disordered metastable solid phase undergoes a 
transition to the crystalline oblique phase as on increas- 
ing temperature system attains sufficient thermal energy 
to overcome local energy barriers and finally undergoes 
a transition to global free energy minimum. 



IV. CONCLUSION 

In this work we have employed computer simulation 
techniques to study the solid-solid collapse transition in 
a model molecular system interacting via both isotropic 
as well as anisotropic interactions. In the case of tem- 
perature induced transition the metastable honeycomb 
solid collapses to high density crystalline oblique phase. 
High density oblique solid has two characteristic length 
scales corresponding to isotropic and anisotropic parts 
of interaction potential. Contrary to the common be- 
lieve and classical nucleation theory that the nuclei of 
new phase are circular (in 2D) and has same order and 
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density as of bulk stable phase, interestingly, we find lin- 
ear strip-like nucleating clusters having significantly dif- 
ferent order and average coordination number than the 
bulk stable oblique phase. Both the order and average 
coordination number of the growing cluster evolve with 
the size and finally converge to the order of the bulk sta- 
ble phase. In the early stage of growth, the cluster grows 
as linear strip followed by branched and ring-like strips. 
Complex geometry of growing cluster is a consequence 
of the delicate balance between two types of interactions 
which enables the dominance of stabilizing energy over 
the destabilizing surface energy of the cluster. The nu- 
cleus of the stable oblique phase is wetted by a layer of 
intermediate order particles which minimizes the surface 
free energy cost for formation of the oblique phase. In the 
case of pressure induced transition the collapsed state is 



a disordered metastable solid. Pathways of pressure and 
temperature induced transitions are significantly differ- 
ent. These observations shed new light on the complex 
nature of nucleation and growth in molecular solid-solid 
transitions as well as on the validity of classical nucle- 
ation theories in complex molecular systems. 
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